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Potassium intercalation in graphite is investigated by first-principles theory. The bonding in 
the potassium-graphite compound is reasonably well accounted for by traditional semilocal density 
functional theory (DFT) calculations. However, to investigate the intercalate formation energy 
from pure potassium atoms and graphite requires use of a description of the graphite interlayer 
binding and thus a consistent account of the nonlocal dispersive interactions. This is included 
seamlessly with ordinary DFT by a van der Waals density functional (vdW-DF) approach [Phys. 
Rev. Lett. 92, 246401 (2004)]. The use of the vdW-DF is found to stabilize the graphite crystal, 
with crystal parameters in fair agreement with experiments. For graphite and potassium-intercalated 
graphite structural parameters such as binding separation, layer binding energy, formation energy, 
and bulk modulus are reported. Also the adsorption and sub-surface potassium absorption energies 
are reported. The vdW-DF description, compared with the traditional semilocal approach, is found 
to weakly soften the elastic response. 



I. INTRODUCTION 

Graphite with its layered structure is easily interca- 
lated by alkali metals (AM) already at room tempera- 
ture. The intercalated compound has two-dimensional 
layers of AM between graphite layers giving rise 

to interesting properties, such as superconductivity.^ 
The formation of an AM-graphite intercalate proceeds 
with adsorption of AM atoms on graphite and absorption 
of AM atoms below the top graphite layer, after which 
further exposure to AM atoms leads the AM intercalate 
compound. 

Recent experiments^' 9 on the structure and elec- 
tronic properties of AM/graphite systems use samples 
of graphite that are prepared by heating SiC crystals 
to temperatures around ~ 1400° CJ£ This heat-induced 
graphitization is of great value for spectroscopic studies 
of graphitic systems, since the resulting graphite overlay- 
ers are of excellent quality^ The nature of the bonding 
between the SiC surfaces and graphite has been explored 
experimentally with photoemission spectroscopy— and 
theoretically 1 ^, with a van der Waals density functional 
(vdW-DF) theory approach that accounts for the van der 
Waals (vdW) forces^^^^ 

Here we investigate with density functional theory 
(DFT) the effects on the graphite structure and the 
energetics and the elastic response when potassium 
is intercalated. The final intercalate compound is 
CsK. The AM intercalate system is interesting in it- 
self and has been the focus of numerous experimen- 
tal investigations i 18 i 19 ' 20 i 21 ' 22 Graphitic systems are also 
ideal test materials in ongoing theory development that 
aims at improving the description of the nonlocal inter- 
layer bonds in sparse systems! 14 ' 23 ? 24 Standard DFT ap- 
proaches are based on local (local density approximation, 
LDA) and semilocal approximations (generalized gradi- 
ent approximation, GGA)^2L2LS for the electron ex- 
change and correlation. Such regular DFT tools do not 
treat correctly the weak vdW binding, e.g., the cohe- 



sion between (adjacent) graphite layers. The failure of 
traditional DFT for graphite makes it impossible to ob- 
tain a meaningful comparison of the energetics in on- 
surface AM adsorption and subsurface AM absorption. 
Conversely, investigations of graphitic systems like CsK 
permit us to test the accuracy of our vdW-DF develop- 
ment work. 

We explore the nature of the bonding of graphite, 
the process leading to intercalation via adsorption and 
absorption of potassium, and the nature of potassium- 
intercalated graphite CsK using a recently developed 
vdW-DF density functional^ This choice of functional is 
essential for a comparison of graphite and CsK properties 
because of the inability of traditional GGA-based DFT 
to describe graphite. We calculate the structure and elas- 
tic response (bulk modulus Bq) of pristine graphite and 
potassium intercalated graphite and we present results 
for the formation energies of the CsK system. 

The intercalation of potassium in graphite is preceded 
by the adsorption of potassium on top of a graphite 
surface and potassium absorption underneath the top 
graphite layer of the surface. In this work we study 
how potassium bonds to graphite in these two parts of 
the process towards intercalation. Our vdW-DF inves- 
tigations of the binding of potassium in or on graphite 
supplements corresponding vdW-DF studies of the bind- 
ing of polycyclic aromatic hydrocarbon dimers, of the 
polyethylene crystal, of benzene dimers, and of poly- 
cyclic aromatic hydrocarbon and phenol molecules on 
graphite^^i^^ 

The outline of the paper is as follows. Section II con- 
tains a short description of the materials of interest here: 
graphite, CsK, and graphite with an adsorbed or ab- 
sorbed K atom layer. The vdW-DF scheme is described 
in Sec. III. Section IV presents our results, Sec. V the 
discussion, and conclusions are drawn in Sec. VI. 
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FIG. 1: (Color online) Simple hexagonal graphite [AA stack- 
ing) and natural hexagonal graphite {AB stacking). The two 
structures differ by that each second carbon layer in AB- 
stacked graphite is shifted, whereas in AA-stacked graphite 
all planes are directly above each other. The experimentally 
obtained in-plane lattice constant and sheet separation of nat- 
ural graphite is (Ref. S3) a = 2.459 A and d c _ c = 3.336 A, 
respectively. 



II. MATERIAL STRUCTURE 

Graphite is a semimetallic solid with strong intra-plane 
bonds and weakly coupled layers. The presence of these 
two types of bonding results in a material with different 
properties along the various crystallographic directions.' 35 
For example, the thermal and electrical conductivity 
along the carbon sheets is two orders of magnitude higher 
than that perpendicular to the sheets. This specific prop- 
erty allows heat to move directionally, which makes it 
possible to control the heat transfer. The relatively weak 
vdW forces between the sheets contribute to another in- 
dustrially important property: graphite is an ideal lubri- 
cant. In addition, the anisotropic properties of graphite 
make the material suitable as a substrate in electronic 
studies of ultrathin metal filmSi 36 ' 37 i 38 ' 39 

The natural structure of graphite is an AB stacking, 
with the graphite layers shifted relative to each other, 
as illustrated in Fig. [TJ The figure also shows hexagonal 
graphite, consisting of AA-stacked graphite layers. The 
in-plane lattice constant a and the layer separation dc-c 
is also illustrated. In natural graphite the primitive unit 
cell is hexagonal, includes four carbon atoms in two lay- 
ers, and has unit cell side lengths a and height c = 2dc-c- 

The physical properties of graphite have been studied 
in a variety of experimenta l 40 ! 41 ' 42 and theoretica l 43 ' 44 
work. Some of the DFT work has been performed in 
LDA, which does not provide a physically meaningful 
account of binding in layered systems ! 15 ' 45 At the same 
time, using GGA is not an option because it does not 
bind the graphite layers. For a good description of the 
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FIG. 2: (Color online) Crystalline structure of CgK show- 
ing the AA-stacking of the carbon layers (small balls) and 
the a/37<5-stacking of the potassium layers (large balls) per- 
pendicular to the graphene sheets. The potassium layers are 
arranged in a p(2 x 2) structure, with the K atoms occupying 
the sites over the hollows of every fourth carbon hexagon. 



graphite structure and nature the vdW interactions must 
be included^ 5 - 

Alkali metals (AM), except Na, easily penetrate the 
gallery of the graphite forming alkali metal graphite in- 
tercalation compounds. These intercalation compounds 
are formed through electron exchange between the inter- 
calated layer and the host carbon layers, resulting in a 
different nature of the interlayer bonding type than that 
of pristine graphite. The intercalate also affects the con- 
ductive properties of graphite, which becomes supercon- 
ductive in the direction parallel to the planes at critical 
temperatures below 1 

The structure of AM graphite intercalation compounds 
is characterized by its stage n, where n is the number of 
graphite sheets located between the AM layers. In this 
work we consider only stage-1 intercalated graphite CsK, 
in which the layers of graphite and potassium alternate 
throughout the crystal. The primitive unit cell of CsK 
is orthorhombic and contains sixteen C atoms and two 
K atoms. In the CsK crystal the K atoms are ordered 
in a p(2 x 2) registry with K-K separation 2a, where a is 
the in-plane lattice constant of graphite. This separation 
of the potassium atoms is about 8% larger than that in 
the natural K bec crystal (based on experimental values). 
The carbon sheet stacking in CsK is of AA type, with the 
K atoms occupying the sites over the hollows of every 
fourth carbon hexagon, each position denoted by a, j3, 
7, or S, and the stacking of the K atoms perpendicular 
to the planes being described by the a/37<5-sequence as 
illustrated in Fig. [2] 



III. COMPUTATIONAL METHODS 

The first-principle total-energy and electronic struc- 
ture calculations are performed within the framework 
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of DFT. The semilocal Perdew-Burke-Ernzerhof (PBE) 
flavor 26 of GGA is chosen for the exchange-correlation 
functional for the traditional self-consistent calculations 
underlying the vdW-DF calculations. For all GGA cal- 
culations we use the open source DFT code Dacapo^ 
which employs Vanderbilt ultrasoft pseudopotentials, 47 
periodic boundary conditions, and a plane- wave basis set. 
An energy cut-off of 500 eV is used for the expansion of 
the wave functions and the Brillouin zone (BZ) of the 
unit cells is sampled according to the Monkhorst-Pack 
schemed The self-consistently determined GGA valence 
electron density n(r) as well as components of the energy 
from these calculations are passed on to the subsequent 
vdW-DF calculation of the total energy. 

For the adsorption and absorption studies a graphite 
surface slab consisting of 4 layers is used, with a surface 
unit cell of side lengths twice those in the graphite bulk 
unit cell (i.e., side lengths 2a). The surface calculations 
are performed with a 4x4x 1 fc-point sampling of the BZ. 

The (pure) graphite bulk GGA calculations are per- 
formed with a 8x8x4 fc-point sampling of the BZ, 
whereas for the CsK bulk structure, in a unit cell at least 
double the size in any direction, 4x4x2 fc-points are used, 
consistent also with the choice of fc-point sampling of the 
surface slabs. 

We choose to describe CsK by using a hexagonal unit 
cell with four formula units, lateral side lengths approxi- 
mately twice those of graphite and with four graphite and 
four K-layers in the direction perpendicular to the layers. 
CsK can also be described by the previously mentioned 
primitive orthorhombic unit cell containing two formula 
units of atoms but we retain the orthorhombic cell for 
ease of description and for simple implementation of nu- 
merically robust vdW-DF calculations. 

In all our studies, except test cases, the Fast Fourier 
Transform (FFT) grids are chosen such that the separa- 
tion of neighboring points is maximum ~0.13A in any 
direction in any calculation. 



A. vdW density function calculations 

In graphite, the carbon layers bind by vdW interac- 
tions only. In the intercalated compound a major part of 
the attraction is ionic, but also here the vdW interactions 
cannot be ignored. In order to include the vdW interac- 
tions systematically in all of our calculations we use the 
vdW-DF of Ref. Uji There, the correlation energy func- 
tional is divided into a local and a nonlocal part, 

E c a £ C LDA + Ef, (1) 

where the local part is approximated in the LDA and the 
nonlocal part E™ 1 is consistently constructed to vanish 
for a homogeneous system. The nonlocal correlation E£ 1 
is calculated from the GGA-based n(r) and its gradients 
by using information about the many-body response of 



the weakly inhomogeneous electron gas: 

Ef = \ [ dr [ dr'n(r)<p(r,r')n(r'). (2) 
2 Jvo Jv 

The nonlocal kernel <f>(r,r') can be tabulated in terms 
of the separation \r — r'\ between the two fragments at 
positions r and r' through the parameters D = (go + 

l'o)\ r ~ r '\/ 2 and & = (Qo- Qo)/(Qo + q'o)- Here q is a 
local parameter that depends on the electron density and 
its gradient at position r. The analytic expression for the 
kernel <fr in terms of D and S can be found in Ref. Il6l . 

For periodic systems, such as bulk graphite, CsK, 
and the graphite surface (with adsorbed or absorbed K- 
atoms), the nonlocal correlation per unit cell is simply 
evaluated from the interaction of the points in the unit 
cell Vq with points everywhere in space (V) in the three 
(for bulk graphite and CsK) or two (for the graphite 
surface) dimensions of periodicity. Thus, the ^-integral 
in Eq. ^ in principle requires a representation of the 
electron density infinitely repeated in space. In prac- 
tice, the nonlocal correlation rapidly converges^! and it 
suffices with repetitions of the unit cell a few times in 
each spatial direction. For graphite bulk the ^-integral 
is converged when we use a V that extends 9 (7) times 
the original unit cell in directions parallel (perpendicular) 
to the sheets. For the potassium investigation a signif- 
icantly larger original unit cell is adopted (see Fig. [2]); 
here a fully converged V corresponds to a cell extending 
five (three) times the original cell in the direction parallel 
(perpendicular) to the sheets for CsK bulk. To describe 
the nonlocal correlations ([2]) for the graphite surface a 
sufficient V extends five times the original unit cell along 
the carbon sheets. 

For the exchange energy E x we follow the choice of 
Ref. [HI of using revPBE^l exchange. Among the func- 
tionals that we have easy access to, the revPBE has 
proved to be the best candidate for minimizing the ten- 
dency of artificial exchange binding in graphite^ 

Using the scheme described above to evaluate E£ , the 
total energy finally reads: 

^vdW-DF _ £.GGA _ ^GGA _|_ ^LDA _|_ jgnl ^ 

where £I GGA is the GGA total energy with the revPBE 
choice for the exchange description and E GGA (£^ DA ) 
the GGA (LDA) correlation energy. As our GGA calcu- 
lations in this specific application of vdW-DF are carried 
out in PBE, not revPBE, we further need to explicitly 
replace the PBE exchange in E GGA by that of revPBE 
for the same electron charge density distribution. 

B. Convergence of the local and nonlocal energy 
variation 

DFT calculations provide physically meaningful results 
for energy differences between total energies To un- 
derstand materials and processes we must compare total 
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energy differences between a system with all constituents 
at relatively close distance and a system of two or more 
fragments at "infinite" separation (the reference system). 
Since the total energy © consists both of a long-range 
term and shorter-ranged GGA and LDA terms it is nat- 
ural to choose different ways to represent the separated 
fragments for these different long- or short-range energy 
terms. 

For the shorter-range energy parts (LDA and GGA 
terms) the reference system is a full system with vacuum 
between the fragments. For LDA and GGA calculations 
it normally suffices to make sure that the charge den- 
sity tails of the fragments do not overlap, but here we 
find that the surface dipoles cause a slower convergence 
with layer separation. We use a system with the layer 
separation between the potassium layer and the nearest 
graphite layer(s) e?c-K = 12 A (8 A) as reference for the 
adsorption (absorption) study. 

The evaluation of the nonlocal correlations E^ 1 requires 
additional care. This is due to technical reasons per- 
taining to numerically stability in basing the eval- 
uation on the FFT grid used to converge the underly- 
ing traditional-DFT calculations. The evaluation of the 
nonlocal correlation energy, Eq. ©, involves a weighted 
double integral of a kernel with a significant short-range 
variation^. The shape of the kernel makes the E£ l eval- 
uation sensitive to the particulars of FFT-type griding, 49 
for example, to the relative position of FFT grid points 
relative to the nuclei position (for a finite grid-point spac- 
ing). 

However, robust evaluation of binding- or cohesive- 
energy contributions by nonlocal correlations can gen- 
erally be secured by a further splitting of energy differ- 
ences into steps that minimize the above-mentioned grid 
sensitivity. The problem of FFT sensitivity of the E™ 1 
evaluation is accentuated because the binding in the E™ 1 
channel arises as a smaller energy difference between siz- 
able E^ 1 contributions of the system and of the fragments. 
Conversely, convergence in vdW-DF calculations of bind- 
ing and cohesive energies can be obtained even at a mod- 
erate FFT grid accuracy (0.13 A used here) by devising a 
calculational scheme that always maintains identical po- 
sition of the nuclei relative to grid points in the combined 
systems as well as in the fragment reference system. 

Thus we obtain a numerically robust evaluation of the 
E^ 1 energy differences by choosing steps for which we 
can explicitly control the FFT griding. For adsorption 
and absorption cases we calculate the reference systems 
as a sum of ^"'-contributions for each fragment and we 
make sure to always position the fragment at the exact 
same position in the system as in the interacting system. 
For bulk systems we choose steps in which we exclusively 
adjust the inter-plane or in-plane lattice constant. Here 
the reference system is then simply defined as a system 
with either double (or in some cases quadruple) lattice 
constant and with a corresponding doubling of the FFT 
griding along the relevant unit-cell vector. 

The cost of full convergence is that, in practice, we of- 



ten do three or more GGA calculations and subsequent 
E^ 1 calculations for each point on the absorption, absorp- 
tion, or formation-energy curve. In addition to the cal- 
culations for the full system we have to do one for each of 
the isolated fragments at identical position in the adsorp- 
tion/absorption cases and one or more for fragments in 
the doubled unit-cell and doubled griding reference. We 
have explicitly tested that using a FFT grid spacing of 
< 0.13 A (but not larger) for such reference calculations 
is sufficient to ensure full convergence in the reported 
(and £; vdw - DF total) energy variation for graphitic 
systems. 

C. Material formation and sorption energies 

The cohesive energy of graphite (G) is the energy gain, 
per carbon atom, of creating graphite at in-plane lattice 
constant a and layer separation c?c-c from isolated (spin- 
polarized) carbon atoms. 

£^G,coh(a, dc-c) = -E-G,tot(a, dc-c) ~ ^G-atom,tot (4) 

where £"G,tot and i?c-atom,tot are total energies per carbon 
atom. The graphite structure is stable at the minimum 
of the cohesive energy, at lattice constants a = ac and 
2dc-c = cg- 

The adsorption (absorption) energy for a p(2 x 2) K- 
layer over (under) the top layer of a graphite surface is 
the difference in total energy [from Eq. for the system 
at hand minus the total energy of the initial system, i.e., 
a clean graphite surface and isolated gas-phase potassium 
atoms. However, due to the above mentioned technical 
issues in using the vdW-DF we calculate the adsorption 
and absorption energy as a sum of (artificial) stages lead- 
ing to the desired system: First the initially isolated, 
spin-polarized potassium atoms are gathered into a free 
floating potassium layer with the structure correspond- 
ing to a full cover of potassium atoms. By this the total 
system gains the energy Ai^K-iayer(aG), with 

A-Ek-laycr (a) = ^K,tot(o) — £"k -atom, tot ■ 

In adsorption the potassium layer is then simply placed 
on top of the four-layer (2 x 2) graphite surface (with 
the K atoms above graphite hollows) at distance dc-K- 
The system thereby gains a further energy contribution 
AEk_q(cIc-k)- This leads to an adsorption energy per 
K-atom 

-Eads(dc-K) = A£ , K -layor(aG) + A^K-G^C-k) • (6) 

In absorption the top graphite layer is peeled off the 
(2 x 2) graphite surface and moved to a distance far 
from the remains of the graphite surface. This process 
costs the system an ("exfoliation") energy — AEq-g = 
-[-Etot.c-G^c-c = c G /2) - £7 tot ,c-G(dc-C -> °°)]. At the 
far distance the isolated graphite layer is moved into AA 
stacking with the surface, at no extra energy cost. Then, 
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the potassium layer is placed midway between the far- 
away graphite layer and the remains of the graphite sur- 
face. Finally the two layers are gradually moved towards 
the surface. At distance 2dc-K between the two topmost 
graphite layers (sandwiching the K-layer) the system has 
further gained an energy A-Ec-k-g(o'c-k)- The absorp- 
tion energy per K-atom is thus 

■Eabs(dC-K) = -A£ C _ G +A£;K-layGr(aG) + A£ C _K-G(dc-K) ■ 

(7) 

Similarly, the CsK intercalate compound is formed 
from graphite by first moving the graphite layers far 
apart accordion-like (and there shift the graphite stack- 
ing from ABA ... to AAA ... at no energy cost), then 
changing the in-plane lattice constant of the isolated 
graphene layers from oq to a, then intercalating K-layers 
(in stacking af3j5) between the graphite layers, and fi- 
nally moving all the K- and graphite layers back like an 
accordion, with in-plane lattice constant a (which has the 
value ac 8 K at equilibrium). 

In practice, a unit cell of four periodically repeated 
graphite layers is used in order to accommodate the 
potassium a/?7i5-stacking. The energy gain of creating 
a (2 x 2) graphene sheet from 8 isolated carbon atoms is 
defined similarly to that of the K-layer: 



Ai?C-laycr (d) — ^C-laycr.tot (a) — 8-Eg- 



atom,tot 



(8) 



The formation energy for the CsK intercalate com- 
pound per K atom or formula unit, inform, is thus found 
from the energy cost of moving four graphite layers 
apart by expanding the (2 x 2) unit cell to large height, 
— Ai?G_ acc , the cost of changing the in-plane lattice con- 
stant from clq to a in each of the four isolated graphene 
layers, 4( Allayer (a) -AI?c-iayer(aG)), the gain of creat- 
ing four K-layers from isolated K-atoms, 4Ai?K-iayer(a), 
plus the gain of bringing four K-layers and four graphite 
layers together in the CsK structure, Ai?c 8 K-acc(a, g?c-k), 
yielding 

Ef oim (a, gJq-k) 



A^G-acc + 4A£ C _ layor (a) - 4A J E C _ 1 ayer(ctG) 



4Ai?K-laycr(a) + A£ Cs K-acc (a, OC-k) 



(9) 



The relevant energies to use for comparing the three 
different mechanisms of including potassium (adsorp- 
tion, absorption and intercalation) are thus -E a ds(^c-K)j 
E&bs(dc-K) and -Ef rm(a, dc-iO at their respective mini- 
mum values. 



IV. RESULTS 

Experimental observations indicate that the intercala- 
tion of potassium into graphite starts with the a&sorption 
of evaporated potassium into an initially clean graphite 



surface. 50 This subsurface absorption is preceded by ini- 
tial, sparse potassium adsorption onto the surface, and 
proceeds with further absorption into deeper graphite 
voids. The general view is that the K atoms enter 
graphite at the graphite step edges^ The amount and 
position of intercalated K atoms is controlled by the tem- 
perature and time of evaporation. 

Below, we first describe the initial clean graphite sys- 
tem, and the energy gain in (artificially) creating free- 
floating K-layers from isolated K-atoms. Then we present 
and discuss our results on potassium adsorption and sub- 
surface absorption, followed by a characterization of bulk 
C 8 K. 

For the adsorption (absorption) system we calculate 
the adsorption (absorption) energy curve, including the 
equilibrium structure. As a demonstration of the need 
for a relatively fine FFT griding in the vdW-DF cal- 
culations we also calculate and compare the absorption 
curve for a more sparse FFT grid. For the bulk sys- 
tems (graphite and CsK) we determine the lattice pa- 
rameters and the bulk modulus. We also calculate the 
formation energy of CsK and the energy needed to peel 
off one graphite layer from the graphite surface and com- 
pare with experiment i^i 



A. Graphite bulk structure 

The present calculations on pure graphite are for the 
natural, AB-stacked graphite (lower panel of Fig.[T]). The 
cohesive energy is calculated at a total of 232 structure 
values (a, dc-c) and the equilibrium structure and bulk 
modulus Bq are then evaluated using the method de- 
scribed in Ref. 

Figure [31 shows a contour plot of the graphite cohesive 
energy variation £g,co1i as a function of the layer sep- 
aration dc-c and the in-plane lattice constant a, calcu- 
lated within the vdW-DF scheme. The contour spacing 
is 5meV per carbon atom, shown relative to the energy 
minimum located at (a, dc-c) — (og, cg/2) = (2.476A, 
3.59 A). These values are summarized in Table Q] together 
with the results obtained from a semilocal PBE cal- 
culation. As expected, and discussed in Ref. 0, the 
semilocal PBE calculation yields unrealistic results for 
the layer separation. The table also presents the cor- 
responding experimental values. Our calculated lattice 
values obtained using vdW-DFT are in good agreement 
with experiment^ - and close to those found from the 
older vdW-DF of Refs.0 and[H, (in which we for as- 
sume translational invariance of n(r) along the graphite 
planes,) at (2.47 A, 3.76 A). 

Consistent with experimental reports^ and our previ- 
ous calculation a 14 ! 15 ' 45 we find graphite to be rather soft, 
indicated by the bulk modulus B value. Since in-plane 
compression is very hard in graphite most of the softness 
suggested by (the isotropic) Bq comes from compression 
perpendicular to the graphite layers, and the value of 
Bq is expected to be almost identical to the C33 elastic 
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FIG. 3: Graphite cohesive energy -Eo.coh {AB '-stacked), based 
on vdW-DF, as a function of the carbon layer separation dc-C 
and the in-plane lattice constant o. The energy contours are 
spaced by 5 meV per carbon atom. 



TABLE I: Optimized structure parameters and elastic 
properties for natural hexagonal graphite (AB-stacking) 
and the potassium-intercalated graphite structure CsK in 
AaApA'yASAa . . . stacking. The table shows the calcu- 
lated optimal values of the in-plane lattice constant a, the 
(graphite-)layer-layer separation d G _c, an< i the bulk modulus 
Bo. In CsK the value if dc-C i s twice the graphite-potassium 
distance dc-K- 

Graphite CsK 
~PBE vdW-DF Exp. PBE vdW-DF Exp 
a (A) 2.473 2.476 2.459 a 2.494 2.494 2.480" 
d C-C (A) > 4 3.59 3.336 a 5.39 5.53 5.35 c 

Bp (GPa) 27 37 de 37 26 A7 de 

"Ref.40. "Ref. 53. c Ref. 4. d Ref. 18. 

e Value presented is for C33; for laterally rigid materials, like 
graphite and CsK, C33 is a good approximation of Bq. 



coefficient j i 18 

We find the energy cost of peeling off a graphite layer 
from the graphite surface (the exfoliation energy) to be 
AE C -g = -435 meV per (2 x 2) unit cell, i.e., -55 
meV per surface carbon atom (Table [Til. A recent 
experiment^ measured the desorption energy of poly- 
cyclic aromatic hydrocarbons (basically flakes of graphite 
sheets) off a graphite surface. From this experiment 
the energy cost of peeling off a graphite layer from the 
graphite surface was deduced to — 52 ± 5meV/atom. 

Our value —55 meV/C-atom is also consistent with 
a separate vdW-DF determination^ of the binding 
(— 47meV per in-plane atom) between two (otherwise) 
isolated graphene sheets. 

For the energies of the absorbate system and of the 
CsK intercalate a few other graphite-related energy con- 
tributions are needed. The energy of collecting C atoms 



to form a graphene sheet at lattice constant a from iso- 
lated (spinpolarized) atoms is given by A.E G _i ayer (a); we 
find that changing the lattice constant a from qq to 
the equilibrium value <zc s k of CsK causes this energy 
to change a mere 30 meV per (2x2) sheet. The contri- 
bution Ai?G-acc is the energy of moving bulk graphite 
layers (in this case four periodically repeated layers) 
far away from each other, by expanding the unit cell 
along the direction perpendicular to the layers. Thus, 
A£ G _ acc = 32A£ G!Coh (a G ,c G /2) - 4A£ c _i aycr (a G ) tak- 
ing the number of atoms and layers per unit cell into 
account. We find the value Ai5 G _ acc = —1600 meV per 
(2 x 2) four-layer unit cell. This corresponds to —50 meV 
per C atom, again consistent with our result for the ex- 
foliation energy, A£" G _ G /8 = —55 meV. 



B. Creating a layer of K-atoms 

The (artificial) step of creating a layer of potassium 
atoms from isolated atoms releases a significant energy 
A-Eic-layer- This energy contains the energy variation 
with in-plane lattice constant and the energy cost of 
changing from a spin-polarized to a spin-balanced elec- 
tron configuration for the isolated atom^ 

The creation of the K-layer provides an energy gain 
which is about half an eV per potassium atom, depending 
on the final lattice constant. With the graphite lattice 
constant a G the energy change, including the spin-change 
cost, is Ai^K-iayer(a G ) = — 476 me V per K atom in vdW- 
DF (— 624 me V when calculated within PBE), whereas 
A-E K -iaycr(ac 8 K) = -473 meV in vdW-DF. 



C. Graphite-on-surface adsorption of potassium 

The potassium atoms are adsorbed on a usual 
ABA . . .-stacked graphite surface. We consider here full 
(one monolayer) coverage, which is one potassium atom 
per (2 x 2) graphite surface unit cell. This orders the 
potassium atoms in a honeycomb structure with lattice 
constant 2a G , and a nearest-neighbor distance within the 
K-layer of a G . 

The unit cell used in the standard DFT calculations 
for adsorption and absorption has a height of 40 A and 
includes a vacuum region sufficiently big that no interac- 
tions (within GGA) can occur between the top graphene 
sheet and the slab bottom in the periodically repeated 
image of the slab. The vacuum region is also large in 
order to guarantee that the separation from any atom to 
the dipole layer— always remains larger than 4 A. 

In the top panel of Fig. 2] we show the adsorption en- 
ergy per potassium atom. The adsorption energy at equi- 
librium is —937 meV per K atom at distance dc-K = 3.02 
A from the graphite surface. 

For comparison we also show the adsorption curve cal- 
culated in a PBE-only traditional DFT calculation. Since 
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FIG. 4: Potassium adsorption and absorption energy at the 
graphite surface as a function of the separation dc-K of the 
K-atom layer and the nearest graphite layer(s) (at in-plane 
lattice constant corresponding to that of the surface, hq). 
Top panel: Adsorption curve based on vdW-DF calculations 
(solid line with black circles) and PBE GGA calculations 
(dashed line). The horizontal lines to the left show the en- 
ergy gain in creating the isolated K layer from isolated atoms, 
A£K-layor( a G)> the asymptote of E ads (d c _ K ) in this plot. 
Bottom panel: Absorption curve based on vdW-DF calcu- 
lations. The asymptote is here the sum A_Ej<;-layer( a G) — 
AEq_q. Inset: Binding energy of the K-layer and the 
top graphite layer ("C-layer") on top of the graphite slab, 
A-Ec-K-G- The dashed curve shows our results when in i?" 1 
ignoring every second FFT grid point (in each direction) of 
the charge density from the underlying GGA calculations, the 
solid curve with black circles shows the result of using every 
available FFT grid point. 



the interaction between the K-layer and the graphite sur- 
face has a short-range component to it, even GGA calcu- 
lations, such as the PBE curve, show significant binding 
(-900 meV/K-atom at d C -K = 2.96 A). This is in con- 
trast to the pure vdW binding between the layers in clean 
graphitei 14 i 15 Note that the asymptote of the PBE curve 
is different from that of the vdW-DF curve, this is due 
to the different energy gains (A-Ex-iayer) in collecting a 
potassium layer from isolated atoms when calculated in 
PBE or in vdW-DF. 

For K-adsorption the vdW-DF and PBE curves agree 
reasonably well, and the use of vdW-DF for this spe- 
cific calculation is not urgently necessary. However, in 
order to compare the adsorption results consistently to 
absorption, intercalation and clean graphite, it is neces- 
sary to include the long-range interactions through vdW- 



DF. As shown for the graphite bulk results above, PBE 
yields quantitatively and qualitatively wrong results for 
the layer separation. 



D. Graphite-subsurface absorption of potassium 

The first subsurface adsorption of K takes place in the 
void under the top-most graphite layer. The surface ab- 
sorption of the first K-layer causes a lateral shift of the 
top graphite sheet, resulting in a A/K/ABAB . . . stack- 
ing of the graphite. We have studied the bonding nature 
of this absorption process by considering a full p(2 x 2)- 
intercalated potassium layer in the subsurface of a four 
layer thick graphite slab. 

Following the receipt of Section IIIII for the absorption 
energy ([7]) the energies AEq-k-g are approximated by 
those from a four-layer intercalated graphite slab with 
the stacking A/K/ABA, and the values are shown in the 
inset of Fig. SJ The absorption energy i? a b s is given by 
the curve in the bottom panel of Fig. [4l and its minimum 
is — 952 me V per K atom at dc-K = 2.90 A. 

To investigate what grid spacing is sufficiently dense 
to obtain converged total-energy values in vdW-DF we 
do additional calculations in the binding distance region 
with a more sparse grid. Specifically, the inset of Fig. 0] 
compares the vdW-DF calculated at full griding with one 
that uses only every other FFT grid point in each direc- 
tion, implying a grid spacing for E^ 1 (but not for the lo- 
cal terms) which is maximum 0.26 A. We note that using 
the full grid yields smaller absolute values of the absorp- 
tion energy. We also notice that the effect is more pro- 
nounced for small separations than for larger distances. 
Thus given resources, the dense FFT grid calculations 
are preferred, but even the less dense FFT grid calcu- 
lations yield reasonably well-converged results. In all 
calculations (except tests of our graphitic systems) we 
use a spacing with maximum 0.13 A between grid points. 
This is a grid spacing for which we have explicitly tested 
convergence of the vdW-DF for graphitic systems given 
the computational strategy described and discussed in 
Sec. III. 



E. Potassium-intercalated graphite 

When potassium atoms penetrate the gallery of the 
graphite, they form planes that are ordered in a p(2 x 2) 
fashion along the planes. The K intercalation causes a 
shift of every second carbon layer resulting into an AA 
stacking of the graphite sheets. The K atoms then simply 
occupy the sites over the hollows of every fourth carbon 
hexagon. The order of the K atoms perpendicular to the 
planes is described by the aj3^5 stacking, illustrated in 
Fig.H 

For the potassium intercalated compound CsK we cal- 
culate in standard DFT using PBE the total energy at 
132 different combinations of the structural parameters a 
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TABLE II: Comparison of the graphite exfoliation energy per 
surface atom, Eq.q/8, graphite layer binding energy per car- 
bon atom, A_Ec_ acc /32, the energy gain per K atom of col- 
lecting K- and graphite-layers at equilibrium to form CsK, 
A-Ec s K-acc/4, and the equilibrium formation energy of CaK, 

Eform ■ 

AE C _ G /S A£ c _ acc /32 AE Cs 

K-acc 

[meV/atom] [meV/atom] [meV/C 8 K] [meV/C 8 K] 

vdW-DF -55 -50 -818 -861 

PBE - - -511 

Exp. -52 ±5" -1236 6 

a Ref. 51. "Ref. 1. 




d c -c [A] 



FIG. 5: Formation energy of CgK, -Ef orm , as a function of 
the carbon-to-carbon layer separation and half the in- 

plane lattice constant, a. The energy contours are spaced by 
20meV per formula unit. 



and dc-c ■ The charge densities and energy terms of these 
calculations are then used as input to vdW-DF. The equi- 
librium structure and elastic properties (B ) both for the 
vdW-DF results and for the PBE results arc then evalu- 
ated with the same method as in the graphite case^ 

Figure [5] shows a contour plot of the CgK formation 
energy, calculated in vdW-DF, as a function of the C-C 
layer separation (cfc-c) and the in-plane periodicity (o) 
of the graphite-layer structure. The contour spacing is 
20 meV per formula unit and are shown relative to the 
energy minimum at (a, c?c-c) = (2.494 A, 5.53 A). 

V. DISCUSSION 

Table [J presents an overview of our structural results 
obtained with the vdW-DF for graphite and CsK. The 
table also contrasts the results with the corresponding 
values calculated with PBE where available. The vdW- 
DF value d C -c = 5.53 A for the C 8 K C-C layer sep- 



aration is 3% larger than the experimentally observed 
value whereas the PBE value corresponds to less than a 
1% expansion. Our vdW-DF result for the CsK bulk 
modulus (26 GPa) is also softer than the PBE result 
(37 GPa) and further away from the experimental esti- 
mates (47 GPa) based on measurements of the C33 elastic 
response^ A small overestimation of atomic separation 
is consistent with the vdW-DF behavior that has been 
documented in a wide range of both finite and extended 
systems.— d5J&LL2L2MiL£l This overestimation results, 
at least in part, from our choice of parametrization of 
the exchange behavior — an aspect that lies beyond 
the present vdW-DF implementation which focuses on 
improving the account of the nonlocal correlations, per 
sc. It is likely that systematic investigations of the ex- 
change effects can further refine the accuracy of vdW-DF 
implementations. 56 In any case, vdW-DF theory calcu- 
lations represent, in contrast to PBE, the only approach 
to obtain a full ab initio characterization of the AM in- 
tercalation process. 

The CsK system is more compact than graphite and 
this explains why PBE alone can here provide a good de- 
scription of the materials structure and at least some ma- 
terials properties, whereas it fails completely for graphite. 
The distance between the graphene sheets upon interca- 
lation of potassium atoms is stretched compared to that 
of pure graphite, but the (K-)layer to (graphite-)layer 
separation, dc-K = dc-c/2 = 2.77 A, is significantly less 
than the layer-layer separation in pure graphite. This in- 
dicates that C 8 K is likely held together, at least in part, 
by shorter-ranged interactions. 

Table Hi] documents that the vdW binding neverthe- 
less plays an important role in the binding and forma- 
tion of CsK. The table summarizes and contrasts our 
vdW-DF and PBE results for graphite exfoliation and 
layer binding energies as well as CsK interlayer binding 
and formation energies. The vdW-DF result for the C 8 K 
formation energy is smaller than experimental measure- 
ments by 31% but it nevertheless represent a physically 
motivated ab initio calculation. In contrast, the C 8 K 
formation energy is simply unavailable in PBE because 
PBE, as indicated, fails to describe the layer binding in 
graphite. Moreover, for the vdW-DF/PBE comparisons 
that we can make — for example, of the CsK layer inter- 
action AEc 8 K-acc — the vdW-DF is found to significantly 
strengthen the bonding compared with PBE. 

It is also interesting to note that the combination of 
shorter-ranged and vdW bonding components in CsK 
yields a layer binding energy that is close to that of 
the graphite case. In spite of the difference in nature 
of interactions, we find almost identical binding energies 
per layer for the case of the exfoliation and accordion in 
graphite and for the accordion in CsK. This observation 
testifies to a perhaps surprising strength of the so-called 
soft-matter vdW interactions. 

In a wider perspective our vdW-DF permits a first 
comparison of the range of AM-graphite systems from 
adsorption over absorption to full intercalation and thus 
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insight on the intercalation progress. Assuming a dense 
2x2 configuration, we find that the energy for potassium 
adsorption and absorption is nearly degenerate with an 
indication that absorption is slightly preferred, consis- 
tent with experimental behavior. We also find that the 
potassium absorption may eventually proceeds towards 
full intercalation thanks to a significant release of forma- 
tion energy. 

VI. CONCLUSIONS 

The potassium intercalation process in graphite has 
been investigated by means of the vdW-DF density func- 
tional method. This method includes the dispersive in- 
teractions needed for a consistent investigation of the 
intercalation process. For clean graphite the vdW-DF 
predicts — contrary to standard semilocal DFT imple- 
mentations — a stabilized bulk system with equilibrium 
crystal parameters in close agreement with experiments. 



Two limits of the absorption process have been inves- 
tigated by the vdW-DF, namely single layer subsurface 
absorption and the fully potassium intercalated stage- 1 
crystal CsK. Here the vdW-DF is shown to enhance the 
(semi-)local type of bonding described by traditional ap- 
proaches. The significant impact on the materials behav- 
ior indicates that the vdW-DF is needed not only for a 
consistent description of sparse matter systems that are 
solely stabilized by dispersion forces, but also for their 
intercalates. 
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